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Abstract. The micropolar fluid mechanics and its transport coeflicients are derived 
from the linearized Boltzmann equation of rotating particles. In the dilute limit, as 
expected, transport coefficients relating to microrotation are not important, but the 
results are useful for the description of coUisional granular flow on an inclined slope. 
(This paper will be published in Traffic and Granular Flows 2001 edited by Y. Sugiyama 
and D. E. Wolf (Springer)). 

1 Introduction 

Micropolar fluids are fluids with micro-structures. They belong to a class of 
fluid with a non-symmetric stress tensor. Micropolar fluids consist of rigid, ran- 
domly oriented (or spherical) particles which have own spins and microrota- 
tions suspended in a viscous medium. The concept of microrotation, originally, 
has been proposed by Cosserat brothers in the theory of elasticity Q. Condiff 
and DahlerQ, and Eringen||^ applied its concept to describe fluids with micro- 
structures in the middle of 60s. Recently, a comprehensive textbook on microp- 
olar fluids has been published^. 

Physical examples of micropolar fluids can be seen in many flelds where all 
of them contain intrinsic polarities. However, we believe that the most inter- 
esting application of micropolar fluid dynamics which would have potential is 
to characterize granular flows In fact, the granular flow is one of flows 
which have micro-structure and rotation of particles. Therefore, it is natural to 
introduce a fluid model which contains an equation for microrotation in addition 
to an equation of velocity. Along this idea, Kanatanij^ formulated micropolar 
fluid mechanics for granular flows, which contains the angle of repose within the 
frame work. Kano et al. |9| have confirmed the quantitative validity of a microp- 
olar fluid model in a chute flow of granular particles from the comparison of 
their simulation of micropolar fluids with their experience. It is worthwhile to 
introduce that the velocity proflle obtained from the micropolar fluid model j^] is 
far from the parabolic curve expected from the conventional Novier-Stokes flow. 

Although we do not know how the micropolar fluid model is relevant in 
other situations of granular materials, it is worthwhile to investigate fundamen- 
tal properties of micropolar fluids from the view of granular physics (see e.g. 
[|o|). Thus, the micropolar fluid mechanics is not a nonsense generalization of 
the Navier-Stokes model, but is a physically relevant model which has many 
applications in physical systems. 
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Recently, Mitarai et al. have analyzed the coUisonal flow of particles on an 
inclined slope. They suggested that the kinetic theory of granular particles is use- 
ful to derive micropolar fluid mechanics, and the fluid equation is quantitatively 
relevant to characterize the granular flow. Their success of the application of 
micropolar fluid mechanics to granular flows makes us understand microscopic 
origin of the micropolar fluid mechanics. Such attempts have been discussed 
within kinetic theory of polyatomic fluids [ ^2^1^ and there are some applica- 
tions to granular fluid [|^, but most of their works are not well accepted in these 
days. 

The purpose of this paper is to summarize our current understanding of 
kinetic theory of micropolar gases. For the sake of simplicity we restrict our in- 
terest to two-dimensional perfectly rough particles in the dilute limit. To avoid 
complication which may be one of reasons for not appealing of the old kinetic the- 
ory of rough particles, we adopt the method of eigenvalue analysis of linearized 
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The organization of this paper is as follows. In section 2 we explain the 
outline of micropolar fluid mechanics and its eigenvalue analysis. In section 3, 
we shortly summarize the classical mechanics of binary collisions of identical 
rough disks. In section 4, we explain the framework of kinetic theory of dilute 
gases of rough disks. We also summarize the relation between the eigenvalues of 
linearized Boltzmann equation and the transport coefficients. In section 5, we 
demonstrate how to obtain transport coefficients in micropolar fluids. In section 
6 we will discuss our results. 



2 Outline of Micropolar Fluid Mechanics and its 
Eigenvalue Analysis 

2.1 Outline of Micropolar Fluid Mechanics 

In this section, we explain the outline of micropolar fluid mechanics. The conven- 
tional fluid mechanics consists of equations of continuity of mass, linear momen- 
tum and energy. In the fluid without microstructure the conservation of angular 
momentum is automatically satisfied but it becomes a nontrivial conservation 
law in the fluid with microstructure. Therefore we need an extra equation of an- 
gular momentum in micropolar fluid mechanics. In coUisional flows we can adopt 
the differential expansion of the strain field in the stress tensor, the couple stress 
and the heat flux. The constitutive equation is similar to that in the Newtonian 
fluid. In this paper we restrict our interest to two dimensional flows. 
The equation of mass conservation is the same as that in usual fluid 

dtp + v- (pu) = (1) 

where dt = d/dt, and p and u = (u^, Uy, 0) are the mass density and the macro- 
scopic velocity field, respectively. 

The continuity equation of linear momentum, on the other hand, becomes 



pDf U = V • (T 



(2) 
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where Dt = 9t + (u- V) is Lagrange's derivative. Here the stress tensor a contains 
the asymmetric part. In coUisional flows in two dimension, the stress tensor Uij 
can be expanded by strains 

= (-_p + CdkUk)Sij + n{diUj + djUi - dkUkSij) 

+firidiUj - djUi - 2ekijWk), (3) 

where Ckij and = wS^^^ are respectively Eddigton's epsilon and the micro- 
rotation (macroscopic spin field). Thus, the micropolar fluid mechanics has the 
spin viscosity fir except for the usual viscosity fj, and the bulk viscosity (. 
The continuity equation of angular velocity is 

nIDtw ^ djCj3 + esijaij, (4) 

where the component of the couple stress for two-dimensional micropolar flow 
is given by C13 = fjLsdxW, C23 = fisdyW and C33 = 0. n and / are respectively 
the number density and the momentum of inertia of each disk. There is the 
relation between n and p as p = nm with the mass of particle m. 
The energy conservation law, in general, can be written as 

pDte = -diQi + aijdiUj + C^jdtWj - ektjCrijWu - F (5) 

where e, F and qi are respectively the energy density, the dissipation rate by 
inelastic collisions and the heat flux. Now we may assume Fourier's law qi = 
—ndiT with the heat conductivity k and the granular temperature. In general, 
the energy density of micropolar fluid consists of two parts, the translational 
energy and the rotational energy. However, the mechanism of the energy transfer 
between two parts is not simple Here, we are only interested in the total 
energy density of two parts. We also assume = to characterize of behavior of 
small dissipations. We will discuss the complete treatment including the energy 
transfer between two parts and the energy loss of inelastic collisions elsewhere. 



2.2 Linearized hydrodynamics and its eigenvalue analysis 

Here let us discuss the behavior of linearized hydrodynamics around the basic 
state p, e and tj = 0: 

p = p + Sp, U = (5U, e — e + Se, (6) 

where U = (u, w) = {uxTUy,'w) is the three dimensional velocity-microrotation 
vector. Assuming the relation between 6e and ST as 6e — CyST with the specific 
heat with constant volume Cy, the set of equations for linearized hydrodynamics 
in Fourier space thus becomes 

dtPq = -ipq ■ u, 

dtUg = -iaq ■ Pq - i/3q • Tg - vq^Uq - (5q(q • u,) + 2w,.q x Wq 

dtWq = -flBQ^Wq + 2flr{iq X Uq - 2Wq) 

dtTq - -Z7q • u, - m^Tq (7) 
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where Wg — (0,0, Wq) is the microrotation vector and the suffix q represents the 
Fourier component of the variable and 
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Mq ■ ^q. 



(8) 

can be summa- 
(9) 



Mq = 






—vq^ —2ivrq 

-2iprq —psq^ — 4/ir 






-if3q 





(10) 



Therefore, the eigen equation can be written as 



(11) 



With the aid of five eigenvectors tp^,, the solution of linearized hydrodynamics 
(^ is represented by 



(12) 



where c^(i) behaves as c^(t) = 0^(0)6"^°*. 

The eigenvalues of eq.(||) is connected with the transport coefficients, where 
the eigenvalues are the solution of 



det{Mq - XII} = 0. 



(13) 



Although the exact solution of (|l^) is difficult to be obtained, the hydrodynamic 
solution near q = can be obtained easily. Following the method introduced in 
standard textbooks (see e.g. |l^) , the five eigenvalues are obtained as follows: 
The first two eigenvalues are 



(14) 



where Cg = {Cp / Cv){dp / d p)t with the specific heat with constant pressure 
Cp. Fs is the rate of sound absorption, which is given hy Fs ^ ^ [C, + p + 



The third mode is 



ri5^ 
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and the fifth mode is 

The fourth mode is the nontrivial mode which represents the relaxation of the 
microrotation as 

A| = p MisH )Q ■ (17) 

We should note that A4 contains a constant which is independent of q. This means 
that the microrotation field cannot be connected with the zero eigenvectors of 
the linearized Boltzmann (or Enskog) equation. 



3 Binary Collisions of Identical Rough Disks 

In this section, we briefiy summarize the result of binary collisions of two- 
dimensional, identical, circular disks obtained by Jenkins and Richman|l7|. 

Here we assume that all variables are non-dimensionalized by the diameter 
of the disk d, the thermal velocity ct = y^Tlj/m, where Tq and m are granular 
temperature and the mass of a particle, respectively. The angular velocity is 
assumed to be non-dimensionalized by lut = \/ 2To // with / — ratP' j%. All 
dimensionless quantities are specified by hat. Thus, we introduce 

a; » c , a; , . 

X = C = , LO = — . (18) 

a Ct ^t 

which means that duj is reduced to 2\f2CjCT- 

Let us consider a binary collision in which the translational velocity Ci, 62 
and the spin cDi, 0)2 priori to a collision to become c'j^, and Qj'i-, ^2- order to 
characterize binary collisions of two identical disks we introduce two coefficients 
of restitution as 

k-v' — —ek ■ V and fc x v' = —fio{k x v) (19) 

where k is the unit vector connecting the center of the particle 1 with that of 2, 
and V is the relative velocity of the points of contact defined by 

V = g + V2z X k{uoi + CJ2) with g = Ci — £2. (20) 

The coefficients of restitution are not constants in actual situation. In particular, 
the tangential restitution coefficient (3q strongly depends on the ratio of the 
tangential velocity to the normal velocity when the Coulomb slip takes place. 
However, they may be regarded as constants in a wide range approximately. The 
coefficient e may take values from to 1, while /3o may vary between -1 and 1. 
The elastic hard core collision is characterized by e = 1 and (3q = —1 but the 
collisions of perfectly rough disks are characterized by e = /3o = 1, where the 
energy is conserved in the binary collisions and there is still the time reversal 
symmetry. 
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In this paper, we are interested in the case e = /3q = I, because this situation 
aUows the simphfied description with time reversal symmetry but the spin can 
play nontrivial roles. In this situation, we have the changes in a binary collision 
of the relative velocity g 



the velocity 



and the spin 



g'-g = -^v- ^fc(g-A:), (21) 



Ci = ~^-^Hk-g), - C2 = I + • g), (22) 



\/2 - ^/2 

uj'i — till — — —{k X v)z, uj'2 ~ UJ2 ~ (k X v)z. (23) 

3 3 

More general results for any e and /Jq can be seen in the paper by Jenkins and 
Richmanfr^. For example, let us present one result for any /3o as 

/2 

Lbl + u)*2 ~u^~U2 = + /3o)^ ■ (k X v), (24) 

OPO 

where w* denotes the spin of i priori to a collision to become a);. Note that there 
is time reversal symmetry for any e and /3o- 



4 Kinetic Theory of Dilute Gases 

In this section, we discuss the kinetic theory of dilute rough disks. Since the 
Boltzmann equation should have modification if there is no time reversal sym- 
metry, here we assume that the disks are perfectly rough, i.e. e = /3o = 1 to 
keep the time reversal symmetry. As will be shown, the result recovers usual 
Navier-Stokes equation in the dilute limit, but in some situation like in flows 
on an inclined rough plate the prediction of the dilute gas kinetics seems to 
be quantitatively useful||ll|. The extension of this analysis for dense particles 
with dissipative collisions will be discussed elsewhere. Note that there are some 
treatments of three dimensional dissipation- less polyatomic fluids [ p"2|p^ based 
on Chapman-Enskog scheme. On the other hand, Jenkins and Richmanp^ de- 
veloped kinetic theory of rough inelastic disks in two dimensions based on the 
Grad expansion. It is difficult to check their argument because their calculation 
is long and complicated. In addition, their method cannot derive micropolar fluid 
mechanics as a closed form. Here we introduce a simpler method of calculation 
which can derive micropolar fluid mechanics and can determine the transport 
coefficients. 

We begin with Boltzmann equation of particles with diameter d: 
dtfi + (ci • V)/i ^djdc2j dLU2 J dkH{g ■ k)g{f[f^ - /1/2} (25) 
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Here we consider the collisions (c*,a;*) — > {c,uj) (c', lj'), but we note {c*,uj*) ~ 
{c',Lo') because of time reversal symmetry ior e — f3Q — 1. In eq.(p5|), we use the 
following notations: k is the common unit normal vector at contact, g = Ci — C2, 
g = |g|, H{x) = 1 for x > and H{x) = for x < 0, = f(r,c[,uj[) and 
/i = /(r,ci,cji). 

Now let us analyze a linearly nonequilibrium situation. The distribution func- 
tion can be expanded as 

/(r, c, Lu, t) = n(r, i)/o(c, + '^'(r, c, uj)) (26) 

where n is the number density, and /□ is the Maxwell-Boltzmann distribution 
function which vanishes in the coUisional integral in eq.(j2^). 

Now we adopt dimensionless quantities as in the previous section. Here we 
only introduce two variables as 

n{x)^n{x)cf, f,^^^M{i) (27) 

where | = (c,6D) and Af(^) = exp(-|^)/7r3/2. 

It is easy to show that the zeroth order equation is given by dn/dx — 0. Thus, 
the distribution function / depends on the spatial coordinate only through ]!'. 

The perturbative equation is now reduced to 

dt^'i + (ci • V)ifi = h j d^i j dkH{g ■ k)gM{i){^!^ +^[-qf^- qf^} (28) 

This equation can be written as 

dt^i + {ci-\I)^i^nL[^i], (29) 

where L\^i\ is the coUisional integral in eq.(p8|). Introducing Fourier transform 
we may rewrite (^ as 

dt^q + iqc^^q = nL[^q]. (30) 
The solution of ( |30| ) can be obtained from the non-Hcrmitian eigenvalue problem 

{L-iqc,)\^'^^>=\'\m> . (31) 



Since we are interested in hydrodynamic behavior of eq.(|3l|), we adopt the 
expansion around q = as 



\1 - \(") ^„2x(2) 



A« ^}}^'> +q}};> +q'Xf> + ■■■. (32) 
Substituting ( ^2| ) into ( ^ ) we obtain 

nLWf>=\f\^f^>. (33) 
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Thus, > is represented by the Hnear combination of five fundamental eigen- 
vectors as: 

5 

I'^i"^ >= E >, (34) 

a' = l 

where \(j)a > satisfies i|(/)Q, >= Xa^\(l)a >, and their expficit expressions are 

104 > = ^, 105 >= - I). (35) 

Since \(j)a > is the degenerated eigenvector, we need to use >. The determi- 
nation of Itf'j"^ > will be discussed later. We also assume that the eigenfunctions 
are orthonormal as 



< 



Therefore we may introduce Aj"'' as 

-xf^n<'F^"^\L\'F^"^>. (37) 
This Xj"' is equivalent to X'"^^ if A^^'' is independent of ^. If A^"' is a function 



of ^, two eigenvalues are different from each other. We believe that A^- plays 



[(0) 

fundamental roles in later discussion. 

With the aid of dSQ) and (pm we obtain the relations at the first order: 



riL ~ Aj 

and 

AW=z<tf'f)|gc,|ip-f >, (39) 

where we use L = Lt and n < w[°^\L = xf ^ < 

From (^ , (|3^ ) and (^) we obtain the second order correction of the eigen- 
value as 

^(2) ^ _ < ^(°)|(jgc, + A«) ^-^{tqc, + xl'^Wf > . (40) 

nL - A^. 

As the eigenvalues A° = are degenerate, we must be careful in starting 
from a proper basis that avoids the appearance of vanishing denominators. As 
in the case of quantum mechanics, we must solve exactly eigenvalue problem 
in the subspace spanned by |0ct >• From the comparison of (^8|) with (34) we 
obtain 

5 

Caa'[< 4'a'\Cx\(f>a" > -KSa,a'] = 0, (41) 
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where A — iX^a'^ /q. The required condition to exist nontrivial solutions of this 
equation is 

det[< <l)i3\c^\(t)c,> -\5c()\=Q (42) 

Thus, we obtain 



Ai = -A2 = iy^, A3 = A4 = A5 = 0. (43) 
For later convenience, we introduce co = Ai. Coining back to eq.(|4]|) we obtain 



1 

71 



> 



> 



I'^r > 



1 

V2 

103 > 

104 > 

2 



h > + 02 > 



^101 > -102 > 



105 > 



^105 > 



-101 > > 



These \^a '' > will be used as the basis of perturbative treatment. 



(44) 



Taking into account (^ we obtain the eigenvalues up to the second order as 
Ai = -icoq + < •ff 



-2 ^ ,T,i")if;^. _.„^J_r^ ._M,r,(o) 



A2 

A3 
A4 

As = < 



icoq + q^ < 'Pi^\{cx + Co)4y{cx + Co)W^"' >, 



1 



,(0) 



1' <€'%^^M'' >, 



Af +9^<ff I 



(0) 



1 

' fiL 

y(0) 



fiL 



(45) 



Basically, the eigenvalues obtained in this section should be equivalent to 
those obtained in section 2 as the linear micropolar hydrodynamics. We have to 
note, however, that the results in this section is written in dimensionless forms 
but the results in section 2 include physical dimensions in terms of the thermal 
velocity ct — \/ 2Tq / m and the diameter of particles d. Thus we obtain the 
relations 



Co 



r. 



m 



'cx-co) — {cx-co)\1'i°\> 
nL 
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TO V TO nL 

2^ -(0)lA ^ A l,r,(") 



-pdCpd— < n">\c,—c,m"> > . (46) 

V TO nL 
5 Evaluation of transport coefficients 

In this section, we evaluate the transport coefficients /x, /i^ and fig among many 
transport coefficients. Actually we can describe incompressible micropolar fluid 
mechanics in terms of these three coefficients. 



5.1 Evaluation of fir 

As was shown in the previous section, to evaluate /i^ we have to obtain A^*^^ . We 
here present the result for any /Jq, since it is possible to obtain A4"'' for any /?o 
in dH). 

Substituting (|4|) into ( p8| ) we obtain 

m - - ^^iy'h c,[c.]+c,m = ^u., (47) 

where 

r-i[u,] ^- J d^C2 J dkHig ■k){g- k)M2ic)z • (fc X g) (48) 

£2^] ^V2Co j d^C2 J dkH{g, • fc)(g • fc)M2(c) (49) 

with M2(c) = cxp[— — Cy]/TT. With the aid of formulae for the modified Bessel 
function Io{z) with zcroth order and the confluent Hypergeometric function 
Fia,b;z): 



and 



Io{z) = - dxe''°^\ f^(3/2,l;c2) = — / drr^h 
Jo Vi" 



(2rc) (50) 



it is possible to show 



C2[u}]^ ^^F{i/2,l-c^)e-^ Lo, and £i[w]=0. (51) 



Therefore, we obtain the eigenvalue 

,(0) _ V^(l + ^o)", 



Ar = - \^ ' e'^ ^(3/2, 1; c^). (52) 
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This result states that the eigenvalue depends on the velocity of particles. 

As mentioned in section 4, if the eigenvalue depends on the velocity, the effec- 
tive eigenvalue Aq defined in ( |37| ) plays important roles. Thus, xf'' is evaluated 
as 

Af = -^^^^ I rf^l5^e-V(3/2,l;c^)^^ (53) 
The result of integration becomes 

Af = ^^^±M^ (54) 



where we use dcce '^'^ -F(3/2, 1; c^) = l/\/2. We note that the result (54) is 
reduced to 



Af = -^n (55) 



in the case of /3o = 1 . 

From (E5|), (Eol) and (p5|) we obtain 



lir = md^J— (56) 

96 V TO 

for e = i3q = 1. Thus, /i^ is proportional to as predicted by Mitarai et 
al.||ll|. Therefore, the effect of microrotation is negligible in usual situations of 
dilute gases. However, as demonstrated by Mitarai et aZ.[0 the shear near the 
boundary produces the relevant situation for microrotation. 



5.2 Evaluation of /i. 

H can be evaluated from the third equation of ( ^ ) . This equation may be rewrit- 
ten as 



J d^c.^Cyx'^yMi^) = - < x^^lr-^ >, (57) 



where \x^'^ >= X^^ \Y^y >= CxCy. ft is defined by 



f,^f^/ipdJ^) (58) 
V m 

and x^y is the solution of 

fiLlx'^y] = c^cy or hL\x''y > . (59) 

Introducing the function j/^*' > vertical to \(j)a >, it is easy to show that 
< f\L\f > is minimum for |/ >= Ix'^v > where |/ >=< > / < 

> 1/^^ >. The minimum principle then becomes 

<x|r-><</|nL|/>=L^l_^^. (60) 
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Here we assume the following expansion: 

m 

ir''>=E/3rV.>, (61) 

where the coefficient /3^^'"'* is determined by the minimum condition of the right 
hand side of (|60|). These conditions are summarized as 

1 ™ 

Y-y = ^ ^(") 6 J ^ < I y-y > = _ ^ ^ ;3(™) Y;^y , (62) 

where MF-^ > and = - < ML|/j >. 

We assume that \ fj > can be represented by 

\fj >^ St'\c')£.cy, (63) 

where S2 (x) is the Sonine polynomial which is defined by 

j:. (-iprq + r + i) 
' ("^-^^r(; + , + i)(r-,)!j!" (^^^ 

with the Gamma function r{x). The Sonine polynomials satisfy the orthonor- 
mality condition dxx'- s\^\x)S\^ \x) = '^^''^''''^^ (^r,r'- In particular, the 
relations sf*^ {x) = 1 is useful for later discussion. From the orthonomality the 
condition of < f^y\(j)a >— J2T f^'i"^^ ^ fj\4'a >= is automatically satisfied. 



■'3 "^J 

Therefore we obtain the expressions for Y^^ and b'^f as 



blf =- I d^MiOsi''-'\c^)c,CyLc,CySt'\c^)- (65) 



It is easy to calculate Yj as 



1 



yr = ihi- (66) 

It is notable that Yj^^ is zero except for Y^^ . 

Since it is known that the expansion in terms of Sonnine's polynomial is fast, 
the result with m — 1 gives a good approximation. Adopting this approximation 
(HI) becomes 

< x'^ir"'' >= -^r"^; Y^y = M^^n- (67) 

Eliminating /J^^-* from ( |67| ) we obtain 

< x^^lr^^ >=-%y- (68) 
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From the definition of the operator L can be written as 

X [c'lxC'ly + c'2Av - ClxCly " C2xC2y\'^ ■ (69) 



where tt — -0 is the angle between g and k. To derive (|69|) we have used that the 
cross section (da /dip with the impact parameter a) is cosTp and the time reversal 
symmetry of collisions. 

From and (^) we get the relation 

From @, (HI) and (||) we obtain 



From comparison of (46) with the definition of /t, the final expression of fi is 
given by 

Note that this result is deviated from the result in ref. . 



5.3 Evaluation of 

The method of calculation of is similar to that of /i. Let us introduce 

Mc = - 1 diCoCxr'^Mii) = - < >, (73) 

where >= a)c^ and \x^^ >= x"^ is the solution of 

nL\x'^^ >= jy^^ > . (74) 
Similar to the previous section fie is represented by 



c2 



f^c = ^ (75) 
in the lowest order approximation. Here Y^^ and are respectively given by 
yr - J d^luJCxfMi^), b't^ diM{i){oCxLuCx. (76) 

It is easy to show 

yr = \- (77) 
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On the other hand, b^'f is rewritten as 



X [c[^L0[ + C2^UJ2 - CixOJl - 02x^2^- (78) 

From (H) and (|2|) we obtain 
Thus becomes 



From and (|7^) we obtain 

3 



Mb = 



d^/^,-X^I^MJ-d^/2^o- (81) 
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It is obvious that the second term is neghgible in the dilute Umit. 

The method of derivation of fis is relatively simple, when we compare the 
method based on Chapman-Enskog scheme in which fi b becomes the correction 
term of higher order ||l2]. 



6 Discussion 



Here, we have demonstrated how micropolar fluid mechanics can be derived 
from the Boltzmann equation. Of course, the result is reduced to Navier-Stokes 
equation in the dilute limit. In this sense, at least, we will have to extend our 
work to the case of Enskog equation if we believe that the concept of micropolar 
fluid mechanics is useful. The effect of dissipation also plays important roles 
though we do not discuss it. 

Roughly speaking the effect of microrotation is localized in the boundary 
layer. Therefore it will be important to analyze simple shear flows. Our prelim- 
inary result suggests that the micropolar fluid mechanics may not be enough to 
discuss the region close to a flat boundary. The success by Mitarai et al. [0 may 
come from their bumpy boundary condition in which there is no Knudsen's layer 
in the system due to random scattering of particles near the boundary. 

There are many problems to be solved. Let us close the paper with our 
perspective whether micropolar fluid mechanics is useful. In the pessimistic view 
the analysis presented here is general nonsense, and micropolar fluid mechanics 
is useless. On the other hand, in the optimistic view, this work will be a milestone 
to discuss the fluid motion with microstructure. Although we are not sure which 
result we will see in the future, we hope that micropolar fluid mechanics is a 
useful concept to characterize the flow of particles. 

The author thanks N. Mitarai and H. Nakanishi for fruitful discussion. 
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